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1. Introduction 

The onset of extragalactic cosmic rays and the end of the Galactic cosmic ray (CR) component 
are subjects of major investigation at this time. However this transition region remains poorly 
understood [p and the conclusions we draw depend rather sensibly on the mass composition and 
on the spectral shape of extragalactic CRs. Of particular importance is the fact that in most models 
of the transition from Galactic to extragalactic CRs one is forced to assume that the flux of the 
extragalactic component is suppressed at sufficiently low energies. Some constraints on the sources 
of ultra high energy CRs (UHECRs) can be obtained as follows [^: for a source with size R 
to be able to accelerate CRs to energy E, the Larmor radius E/ZeB must be smaller than the 
size R (here Ze is the charge of CRs and B is the magnetic field). This implies that B > E jZeR. 
Since the source must be able to guarantee at least a magnetic energy flux vB^/An, where v is 
the velocity of the accelerator, one easily gets a lower limit on the source luminosity: L>3.2x 
) j8 erg/s where jS = v/c is the dimensionless velocity in units of the speed of 
light. The required luminosity becomes even higher in the case of relativistic sources [^. The 
general conclusion is that the sources of UHECRs have to be very luminous. 

Being interested in investigating the escape of CRs from their sources, here we will not con¬ 
sider the problem of acceleration any longer and we will concentrate on the propagation outside 
the host galaxy inside which the actual source is located. On the other hand, for the sake of brevity 
we will refer to such galaxy as the source. 

Eet us assume, for simplicity, that CRs leave their host galaxy with an injection spectrum 
q{E) oc £^2 yp jQ some maximum energy Emax- The differential number density of CRs streaming 
out of such sources can then be written as: 


ncR{E,r) 


q{E) 

ATlr^c 


Lcr E-^ 
A ATir^c 


^l.lx 10-'%4 ^Mpc cm-^GeV-i , (1.1) 

where we have adopted A = ln(Emax/Emin) ~ 25 and Lcr = 10'*^E44 erg/s, energies are in GeV and 
distances in Mpc. The source is assumed to be located in an average position of the intergalactic 
medium (IGM) with a density of baryonic gas is nt, = Q-bPcr/nip = 2.5 x 10^^ (u^) cm^^. We 
also assume that there is a cosmological magnetic field with a strength Bq = G and a 

correlation scale of the order of Ag ~ 10 Mpc, so that on scales smaller than Ar, the field can be 
considered as oriented along a given £ direction. Hence the Alfven speed is va = Eq/ \/^^^hPcr ~ 
44 B_i3cm s^^ 

Assuming that CRs escaping the host galaxy are all protons (positively charged particles), one 
can easily see how this situation leads to the formation of an electric current directed outwards. 
In order to avoid charge accumulation anywhere in the Universe, one is forced to require a return 
current that exactly compensates the CR current. This situation, similar to that the takes place 
near supernova remnants, is known to be able to potentially lead to instabilities. Investigating such 
instabilities and their implications is the main purpose of this paper. 
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2. Calculations 


The current of CRs streaming away from their source into the IGM can be written as 


7cr = encR{> E)c = 


eLcRE 

AnAr'^ 


(2.1) 


where E is the minimum energy of the particles constituting the current. The validity of this 
expression is restricted to the case of vanishing background magnetic field, or a case in which the 
Larmor radius of the particles is r/, 3> Ag, but we shall see that the above estimate for the current 
density turns out to hold also in the diffusive regime. 

A current propagating in a plasma can give rise to instabilities of different types. Granted that 
the current carrying particles are well magnetised (va > JcR/{^n},), which from Eq. ^ is seen to 


imply Bq > 2 X 10 


^GeV''Mpc 


) G, the fastest growing instability arises when the condition 


JcrE > 


ccBq 

An 


( 2 . 2 ) 


is satisfied. This condition, which is the standard one for the development of non-resonant modes 
of the streaming instability [|1], is easily seen to be equivalent to the requirement that the energy 
density of CRs be locally larger than the energy density in the form of pre-existing magnetic field, 
B^/An. For q{E) oc E^'^, fhis requiremenf becomes independenf of E and, using Eq. |2.1[ , can be 
simply formulafed as: 


. 1/2 


r < rinst = 3.7 X 10'’ 


B 


13 


Mpc. 


(2.3) 


When Eq. is salisfied, fhe fasfesl growing modes have a wavenumber k^ax that reflects the 
equilibrium between magnetic tension and 7 cr x 8B force on the plasma, namely kmax^o = ^^cr, 
and their growth rate is: 


Imax — ^max^A — 


/ 471 7cr 

Ufo/Mp C 


(2.4) 


independent of the initial value of the local magnetic field Bq- The scale of fhe fasfesf growing 
modes kj^ax much smaller fhan fhe Earmor radius of fhe particles dominafing fhe currenf (fhis is 


enfailed in fhe condition for growfh, Eq. |2.2| ), fherefore fhey have no direcf influence on parficle 
scattering. This conclusion is however changed by fhe non-linear evolution of fhe modes. As long 
as fhe insfabilify develops on small scales, if cannof affecf fhe currenf, hence one could freaf fhe 
fwo as evolving separately. A fluid elemenf will be subjecf fo a force fhaf is basically ~ Jcr8B/c'. 
ifs equation of motion is p{dv/dt) ~ 2/cr 55> wifh 5B{t) = 5Boexp (ymaxO- estimate 

fhe velocify of a fluid elemenf as v ~ (5B(t)/cR)/(cp7max)> which upon integration in time leads 
fo an estimate of the mean fluid displacemenf as Ar ~ (^B(t)7cR)(cPT ^ax )- The safurafion of fhe 
insfabilify is achieved when fhe displacemenf equals fhe Earmor radius of parficles in fhe currenf, 
as calculated in fhe amplified magnetic field, E / e8B\ when fhis condition is fulfilled, scattering 
becomes efficienl and fhe currenf is desfroyed. This simple criterion refurns fhe condifion: 


5b2 JcrE 


An 


ce 


= ncR{>E)E. 


(2.5) 
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Since ncR(> E) oc in the case considered here, the saturation value of the magnetic field is 
independent of the energy of particles in the cuiTent driving the instability. A somewhat lower 
value of the saturation can be inferred following the criteria discussed in [Q], where the non-linear 
increase of the wavelength of the fastest growing modes is followed. Using such a prescription 
our saturation magnetic field would be ~ 10 fimes smaller. Eq. ^ expresses fhe condifion of 
equiparfifion befween fhe CR energy densify and fhe amplified magnetic pressure, a condifion fhaf 
is offen assumed in fhe liferafure wifhouf jusfificafion, and fhaf here arises as a resulf of fhe physics 
of magnetic amplification ifself. 

The field sfrengfh, as a function of fhe disfance r can now be written as 


5B{r) = 3.7 X 10 Gauss . 


( 2 . 6 ) 


This rafher sfrong magnefic field will develop over disfances r from fhe source fhaf safisfy Eq. 
and under fhe addifional condifion fhaf fhe growfh is fasf enough so as fo reach safurafion in a frac- 
fion of fhe age of fhe universe, to (numerical simulafions of fhe insfabilify show fhaf safurafion 
occurs when Jmax'^ ~ 5 — 10). This latter condifion reads /maxfo ^ 5 and franslafes info: 


If fhe condifions expressed by Eqs. 2.3 and 2.7 are fulfilled, fhen fhe magnefic field can be esfimafed 
as in Eq. ^ and since dB » Bq and fhere is roughly equal power af all scales, if is reasonable fo 
assume fhaf parficle propagafion can be described as Bohm diffusion in fhe magnefic field 5B, so 
fhaf 


D{E,r) = 9 X lO^'^EoeV ?'Mpc E 44 ^^cm 


( 2 . 8 ) 


The inifial assumpfion of ballisfic propagafion of CRs escaping a source leads fo conclude fhaf 
particles would produce enough furbulence fo make fheir mofion diffusive. The diffusion fime over 
a disfance r from fhe source can be estimated as Xd{E, r) = /D{E,r) 3.3 x lO'^rMpc ^cev 

from which follows that particles can be confined within a distance r from the source for a time 
exceeding the age of the Universe, if their energy satisfies the condition: 


E < Econf = 2.6 X lO'^ rMpc GeV . 


(2.9) 


One might argue that this conclusion contradicts the assumptions of our problem: for instance 
the density of particles in the diffusive regime is no longer as given in Eq. O- This is certainly 
true, but the current that is responsible for the excitation of the magnetic fluctuations remains the 
same, as can easily be demonstrated: for particles with energy E > Econf in Eq. ^ , and assuming 
that energy losses are negligible, quasi-stationary diffusion can be described by the equation 


li. 

dr 


r^D{E,r] 


dn 


q{E) 

Anr'^ 


5{r) 


( 2 . 10 ) 


where q{E) is the injection rate of particles with energy E at r = 0. This equation is easily solved 
to provide the density of CRs at distance r from the source: 


n{E,r) 


9 

?,KrD{E,r) 


( 2 . 11 ) 
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Clearly, by definition of diffusive regime, the density of particles returned by Eq. 2.11| is larger than 
the density in the ballistic regime, Eq. O- However, the current in the diffusive regime is 


dn q{> E) 




(2.12) 


which is exactly the same current that we used in the case of ballistic CR propagation (see Eq. [2.1[ ). 


This is a very important and general result: the magnetic field in Eq. is reached outside a CR 
source independent of the mode of propagation of CRs, since it is only determined by the current 
and not by the absolute value of the CR density. Clearly the particles that are confined within a 
distance r around the source do not contribute to the CR current at larger distances. 


3. Results and implications 

The energy below which CRs can be considered as confined in the source vicinity can be 
calculated by taking into account all the three conditions that need to be imposed to guarantee such 
a confinement, namely Eq. (existence of fastly growing modes), Eq. (growth rate faster than 
the expansion of the universe) and Eq. (confinement). The first condition yields a limit on the 
distance from the source that is easy to satisfy, unless the strength of the background magnetic field 
is increased by several orders of magnitude, in which case however other complications arise (see 
discussion below). 

The other two conditions lead to the constrain 


Ecut~10'GeVL: 


2/3 

'44 • 


These particles are confined within a distance from the source: 

r-conf ~ 3.8 Mpc 


(3.1) 


(3.2) 


_ _ :9.6x 

G. It is noteworthy that both the size of the confinement region and the magnetic field 


Within such a distance the magnetic field is as given by Eq. 2.6 and larger than 5B(rconf) 

- 1/3 

-44 


depend weakly upon the CR luminosity of the source, respectively as and . Hence we can 
conclude that magnetic fields at the level of 0.1 — InG must be present in regions of a few Mpc 
around the sources of UHECRs. As a consequence, the spectrum of CRs leaving these sources and 
eventually reaching the Earth must have a low energy cutoff at an energy Ecut- This kind of cutoff 
has been often postulated in the literature in order to avoid some phenomenological complications 
that affect models for the origin of UHECRs. Eor instance, a low energy cutoff is required in the 
dip model [||, ^ to describe appropriately the transition from Galactic to extragalactic CRs. This 
feature is usually justified by invoking some sort of magnetic horizon in the case that propagation of 
UHECRs is diffusive in the lower energy part of the spectrum [^. A similar low energy suppression 
of the CR flux is required by models with a mixed composition [||]. In the calculations illustrated 
above, the presence of nuclei is readily accounted for, provided the current is still produced by 
protons (assumed to be the most abundant specie). In this case, the value of Ecut is simply shifted to 
Z times higher energy for a nucleus of charge Z. The model proposed here provides a mechanism of 
self-confinement of CRs close to their sources that does not need artificial assumptions on diffusive 
propagation. 
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As discussed above, the condition that guarantees the existence of non-resonant modes (Eq. 2.3) 


is easily satisfied, unless the background magnetic field reaches Bq 9.6 x 10^^^ . However, 

when fhis happens fhe calculafions above break down for anofher reason: CRs can freely sfream 
from fhe source only if fheir Larmor radius in fhe pre-exisfing magnefic field is much larger fhan 
fhe assumed coherence scale of fhe field, namely 


lO^GeV 


13 ^ 10 - 


(3.3) 


When Eq. ^ is nol fulfilled, namely when fhe background field is relatively sfrong, fhen fhe 
propagafion of CRs from fhe source becomes infrinsically one dimensional, which implies fhaf fhe 
densify of particles should be written as 

QiE)t 2Q{E) 


ncK{E,r) = 


(3.4) 


Tlr^vt kRc{E)^c^ 

where we used fhe fad fhaf fhe mean velocify of fhe particles is v = c/2 for a disfribufion of parficles 
fhaf is isotropic on a half plane and we assumed that particles spread in the direction perpendicular 
to the background field by a disfance equal fo Rc{E) = ma.x{rL{E),Rs) wifh R^ fhe source size and 
ri{E) fhe Earmor radius of particles of given energy E. Eor a given source size > Rs as soon as 
^GeV >9X10^5 _io(Ri/100k/rc). Af energies larger fhan fhis: 

ncRiE,r) ^ 47L44EGeVdocm-3GeV-\ (3.5) 

If pai'ticles with energy > £" are able to reach a given location, the current at that location is 


c , , EQ(E) 

JcR ^e-EncR{E,r) =e- 


Kr, 


which resulfs in non-resonanf growfh of fhe field for 


= 3.5 X lO^GeVL, 


1/2 
'44 ’ 


independenf of Bq, and in a safurafion magnefic field 

5B ^ 0.7 B 


10 


Gauss. 


(3.6) 


(3.7) 


(3.8) 


This value of fhe magnefic field is apparenfly very large and reflecls fhe very large densify of 
particles af low energies in fhe proximity of fhe source, as due fo fhe reduced dimensionalify of fhe 
problem. However one should nofice fhaf fhe value is normalized fo fhe densify of GeV particles, 
which only live in fhe immediafe vicinity of fhe source and generafe small scale fields fo which 
high energy parficles are almosf insensitive. Af Mpc scales, where only high energy parficles can 
reach, fhe field sfrengfh is much lower as we discuss below. 

Assuming again fhaf fhe diffusion coefficienl is Bohm-like, one can wrife: 

D{E,r) = 4.8 X 10^'’ E^ev ^^44 cmVs, (3.9) 

which leads fo an esfimafe of fhe diffusion time: Tdiff = 6.2 x 10^^ ^cev ^ 44 ^ ®-io ''Mpc y*"- Requir¬ 
ing fhaf parficles reach fhe location af disfance r in a fime shorfer fhan fhe age of fhe universe, we 
fhen obfain: 

Aonf « 0.5 BZf Mpc. (3.10) 
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Following the usual procedure, one can calculate the growth rate of the fastest modes: 

l4n ecncR{> E) 


Ymax — 


= 1.9x 10'%4fiiio£'( 


s-' 

10-^GeV ^ ’ 


(3.11) 


Pb 2 

and impose the condition that /maxlo > 5, which reads: 

(3.12) 

The intersection of all the conditions listed above leads to conclude that particles with energies 

(3.13) 


E < ^growth ~ 5.3 X 10“ GeVL^f B- 


E < Bcut = 2.2 X lO^GeV Ajo 


will be confined within a radius 

^■conf ~ 10 Mpc Aio . (3.14) 

The amplified magnetic field af such disfance is 

5B^3x 10^‘^G• (3-15) 

We emphasise again fhaf fhe resulfs illusfrafed bofh in fhe case of 3-d (lower Bq) or 1-d propa¬ 
gation (higher Bq) are only sensitive fo fhe CR currenf, and hence insensitive fo whefher parficle 
propagafion is ballistic or diffusive. 


4. Summary 


The sfrong effecl fhaf CRs have on fhe environmenf in which fhey propagafe has long been 
invesfigafed in fhe confexf of particle accelerafion af shock waves IH] and is accompanied by 
observafional consequences []T^ [^, such as spafially fhin rims of enhanced X-ray synchrofron 


emission (see [ ]12| ] for a review). Here we invesfigafed fhese processes in fhe immediafe proxim- 
ify of fhe sources (or hosf galaxies of fhe sources) of UHECRs, where fhe CR currenf generafes 
insfabilifies fhaf can change of way CRs propagafe. 

For weak values of fhe sfrengfh of fhe pre-exisfing magnefic field Bq (say < 10^'° G), in fhe 
absence of non-linear phenomena, CRs propagate in approximately straight lines. The resulting 
electric current leads to the development of a Bell-like instability, that modifies fhe propagafion of 
particles fo be diffusive: we find fhaf particles wifh energy < 10^GeV are confined inside a 
disfance of Mpc from fhe source for fimes exceeding fhe age of fhe Universe, fhereby 

infroducing a low-energy cutoff af such energy in fhe specfrum of CRs reaching fhe Earth. Since fhe 
confinemenl disfance is weakly dependenf on fhe source luminosity, we conclude fhaf a region wifh 
~ nG fields should be presenf around any sufficienlly powerful CR source. If larger background 
magnetic fields are presenf around fhe source, fhe gyration radius of fhe particles can be smaller 
fhan fhe coherence scale of fhe field, and in Ibis case CR propagation develops in basically one 
spafial dimension. Eor a coherence scale of 10 Mpc, CRs are confined in fhe source proximity for 
energies E <2x B^J^^ Aio GeV. 


Whefher nafure behaves in one or fhe ofher way depends on fhe poorly known value of Bq. 


Earaday rofafion measures [13] provide weak, model dependenf upper limifs fo Bq in fhe nG range. 
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while lower limits can be imposed based on gamma ray observations of distant TeV sources [|^, 

{Bo > 10-1’ G). 

The physical prescription adopted here leads to estimating the strength of the self-generated 
magnetic field 5B in the source proximity at the level of equipaitition with the energy in the form 
of escaping cosmic rays, independent of the value of the pre-existing field, Bq. A weak dependence 
on Bq was insfead found for fhe safurafion level in [Q], which in our case would lead fo SB abouf 
~ 10 times smaller for small values of Bq, fhereby reducing fhe energy below which CRs are 
confined in fhe source proximify. Undersfanding fhe dynamics of fhe magnefic field amplificalion 
and safurafion is clearly very imporfanf. One could fesf fhe amplificalion mechanism in fhe case 
of supernova remnanl shocks: in Ihis case fhe safurafion criterion used here franslales lo SB 
s/AhwckVs/c, wilh Vs the velocity of the supernova blast wave and wcr is the energy density in 
accelerated particles. Applying this criterion, we obtain an estimate of the magnetic field which is 
in good agreemenf wifh fhaf measured in young galacfic SNRs On fhe ofher hand, due fo fhe 
relalively small value of SB/Bq, fhe safurafion provided by [Q] would refurn a value of SB only a 
factor ~ 2 smaller, loo small a difference to discriminate belween fhe Iwo eslimales. The lesling is 
Ihen leff lo numerical experimenfs sludying fhe propagation of a currenl of energetic particles in a 
low densify, low magnefic field plasma: hybrid simulations wilh Ihis aim are currenlly ongoing. 

The phenomenon of CR confinemenl illuslraled here has profound implicafions for fhe de¬ 
scription of fhe Iransifion region belween Galacfic and exlra-galaclic CRs [Q, ||]. If is ralher 

lanlalising lhal fhe culoff obfained here as due fo self-lrapping is in fhe same range of values fhaf 
have previously been invoked in fhe lileralure based upon phenomenological consideralions. 
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